Monte Carlo Test of the Classical Theory for Heterogeneous Nucleation Barriers 
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Flat walls facilitate the condensation of a supersaturated vapor: Classical theory of heterogeneous 
nucleation predicts that the free energy barrier AF£ et which needs to be overcome for the formation 
of sphere-cap shaped nucleation seeds is smaller than the barrier AF£ om for spherical droplets in 
the bulk by a factor < f(6) < 1, which only depends on the contact angle 9. In this letter 
we compute both AFj* om and AFj* e t from Monte Carlo simulations and test the theory for the 
lattice gas model (for which 8 can be readily controlled). Even though the theory is only based on 
macroscopic arguments, it is shown to hold for experimentally relevant nanoscopic nucleation seeds 
(20 < AF^ OIa /kBT < 200) if (independently estimated) line tension effects are considered. 

PACS numbers: 64.60.an, 64.60Q-,64.70.F-,68.08Bc 
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Nucleation [lj-Q is a ubiquitous process which gov- 
erns such diverse phenomena as the formation of rain 
drops or snow flakes in the atmosphere, the crystalliza- 
tion of proteins and even the formation of industrially 
relevant polymer foams. Even though considerable em- 
pirical knowledge has been gathered, e.g., on which ma- 
terials are suited as seed particles for making clouds rain 
Q or as nucleation agents to start precipitation in met- 
allurgy [7] , many questions remain. Reasons for this par- 
tial shortcoming are manifold: Experiments are very dif- 
ficult to control. Homogeneous nucleation in the bulk 
needs to be separated from heterogeneous nucleation at 
walls or impurities, and small changes in parameters typ- 
ically lead to large changes in corresponding nucleation 
rates. Nucleation sites often only consist of a few hundred 
or thousand atoms and can only be observed indirectly. 
Unfortunately, theoretical progress is also hampered by 
these difficulties. Simulations on the other hand offer full 
control and are at least in principle able to bridge the gap 
between theory and experiment and, indeed, sig nificant 
progress has been achieved in recent years [§4l(J. 

According to Classical Nucleation Theory [if, the 
formation of a (spherical) nucleation seed in the bulk is 
understood in terms of two competing factors: a volume 
term which seeks to expand the seed and an opposing 
surface term F s — 4irR 2 j: 



AF(R) = -Afi(pi 
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R is the radius of the droplet (or bubble), 7 the in- 
terfacial tension (of a macroscopically flat interface), 
A/i = /i — [i coex the difference in chemical potential rela- 
tive to the coexistence value, and pi <v the densities of the 
coexisting phases. The free energy barrier AF* can now 
be obtained by determining the maximum of Eq. ([1]) at 
R* . Interestingly, at R* we also obtain 



AF* = ^F S (R* 



(2) 



from the knowledge of the surface free energy of a droplet 
or bubble. Note that the spontaneous formation of a crit- 
ical nucleus of the new phase in the bulk requires a free 
energy barrier AF£ om of the order 20 < AF£ om /k B T < 
200. Compared to thermal fluctuations (that one can 
easily detect e.g. by scattering experiments) this is a 
rare event. 

The presence of a wall facilitates nucleation and re- 
duces the barrier. For macroscopic droplets (see Fig. [TJ, 
the contact angle 9 between the droplet and the wall is 
given by the competition of the vapor-liquid (7^), wall- 
vapor ( 7 iuu) and wall- liquid (j w e) interface free energies 
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j v e cos 6 = ~/ wv - 7^, 



(3) 



for incomplete wetting conditions (j v i > j wv —j w e) 
Eq. §3§ is known as Young's equation and can be de- 
rived by considering the forces acting at the three phase 
boundary. For small droplets, a correction due to the line 
tension (r) of the contact line arises [lij (which has the 
length 27rr, where r = Rsm8). 

Considering Eq. Q and the geometry of the sphere- 
cap (see caption of Fig. [1]) we obtain Turnbull's estimate 
for the surface free energy [l3l - [l5j : 
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f(9) = (1 -cos0n2 + cos6>)/4. 



(4) 



Intriguingly, f(0) = V sphere ^ cap /V d ropiet- Therefore, 
Eq. @ also holds for the heterogeneous case and we get: 
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Hence, the nucleation barrier can directly be inferred 



The nucleation barrier in the heterogeneous case is sim- 
ply reduced by a factor of f(9). Note that this derivation 
is solely based on macroscopic considerations and has to 
our knowledge never actually been tested. In this let- 
ter we will demonstrate with Monte Carlo simulations 
of an Ising lattice gas model that this formula can be 
applied to nanoscopic systems if line tension effects are 
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FIG. 1: a) Schematic drawing of the system which we use 
to study stable wall-attached droplets in thermal equilibrium. 
The sphere-cap shaped droplet has height h and covers a circle 
of radius r at the wall, with r — RsinO, h = R(l — cos 9). R 
denotes the radius of curvature and 9 the contact angle. For 
R < oo, the density p' v of the vapor which coexists with the 
liquid droplet (at the same chemical potential p, > /j cocx ) ex- 
ceeds p v (Gibbs-Thomson effect [H,[lJ])- b) Typical snapshot 
picture of a lattice gas system at ksT/J = 3.0, p = 0.065, 
L = D = 40 (measured in units of the lattice spacing), and 
Hi/ J = (9 = 90°). Occupied lattice sites are highlighted 
by dots. 



considered. This provides a stringent test for both the 
classical theory of homogeneous and heterogeneous nu- 
cleation. For this purpose we propose a new method 
Hill to measure surface free energies (and hence nu- 
cleation barriers) of liquid droplets in the bulk and of 
droplets attached to walls over the experimentally rele- 
vant range. No "atomistic" identification of which par- 
ticles belong to the droplet and which to the vapor is 
required. No bias potential to stabilize droplets of a par- 
ticular size 0] is needed, and fluctuations of the droplet 
are not constrained. We independently determine the 
droplet volume and the chemical potential that charac- 
terizes the (stable) equilibrium between the small droplet 
and surrounding vapor in the finite system and derive an 
expression for the line tension. 

We test TurnbuH's formula © with Monte Carlo 
simulations of the standard nearest neighbor (Ising) 
lattice gas model on the simple cubic lattice: H = 



—JJ2SiSj,Si.j = ±1. Phase coexistence between sat- 
urated vapor at density p v and liquid pi occurs at a 
(known) chemical potential /i C ocx- (In magnetic nota- 
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FIG. 2: Contact angle 9 as derived from Eq. ^ plotted vs. 
Hi/ J at ksT / J = 3.0 for two lattice sizes (full and broken 
curves). Symbols show the prediction if the estimates for t 
from our analysis of the droplet excess free energies and r — 5 
are considered (see below). 



tion, the magnetic field H = (// — /i cocx )/2, and pi, p v 
are related to the spontaneous magnetization m cocx as 
p v = (1 - m cocx )/2, pi = (1 + TO coox )/2.) Bulk sys- 
tems are studied in cubic simulation boxes with periodic 
boundary conditions in all directions. For the heteroge- 
neous case we choose a L x L x D geometry with pe- 
riodic boundary conditions in x and y direction. Sur- 
face fields H\ and Hp = —Hi act on (the first layer of) 
the two free L x L surfaces. The temperature is set to 
ksT/J = 3.0 to stay away from the bulk critical tem- 
perature (at ksT/J w 4.51) and above the roughening 
transition temperature [18| (at ksT/ J w 2.45). Hence, 
the correlation length £ is less than a lattice spacing and 
the interface free energy does (to a good approximation) 
not depend on the interface orientation. 

Wetting for the Ising model has been thoroughly 
studied before [19| and by varying the surface field 
H\ we can control the contact angle (Fig. [5]). Note 
that the local surface layer magnetization mi = 
— (df s (T,H,Hi)/dHi)T,H [20] (f s is the surface excess 
free energy per spin). With f s (T,H — )• 0~,Hi) — j wv , 
f s (T,H — ¥ + ,Hi) — j W £ and Young's equation ([3]), we 
obtain 9 from thermodynamic integration 



Hi 



cos 9 = (l/7„<) / (m D - mi)dH[ . 



(6) 



Here, we exploit the Ising symmetry f s (T,H — > 
0~,Hi) = f s (T,H -> 0+,-iJi) to use the local mag- 
netization rriD at the other surface at which Hjj = —Hi 
acts. "f v i = 0.434 was given by the very accurate esti- 
mates of Hasenbusch and Pinn [2l[ . 

In the following we describe how to obtain the surface 
free energy as a function of the droplet radius. First, 
we measure the chemical potential (in the canonical en- 
semble) as a function of density (Fig. by applying an 
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FIG. 3: Chemical potential as a function of density. The den- 
sity of the sphere-cap shaped droplet at overall density p is 
given by p[, the density of the environment by p' v . AF S is ob- 
tained by integrating the chemical potential over the density 
from p' v to p. 



adaptation |16|, [17| of the Widom test particle insertion 
method (2^. In the flat region of the isotherm depicted 
in Fig. [31 a single droplet is present in the simulation box, 
which will be used for the determination of the surface 
free energy. At small densities the ascending branch of p 
represents pure vapor (with no single large droplet being 
present), and the subsequent rapid decrease signifies the 
system size dependent transition towards a single droplet 
phase f23T - |26j ] . At larger densities, the droplet region is 
delimited by a transition towards a cylindrical droplet 
which is stabilized by the boundary conditions. Note 
that for each system size, only a certain range of droplet 
sizes can be stabilized and we have to simulate a whole 
range of box sizes to cover all radii. 

As indicated in Fig. [3] for each droplet state at den- 
sity p, there exists one vapor state at p' v and one liquid 
state at p[ which have the same chemical potential. Equi- 
librium conditions require that the chemical potential is 
constant throughout the box. Therefore, the droplet has 
density p\ and is surrounded by an environment with den- 
sity p' v . The surface free energy of the droplet is simply 
given by the difference in free energy between the droplet 
and the vapor state and can be determined by thermo- 
dynamic integration: 



(7) 



F s (p) = / ii(p')dp'- 

•'P'v 

The number of particles in the box is conserved: 



Pv(V ~ Vdrop) + p'lVdrop = PV. 



(8) 



If we further require that the (average) droplet is spher- 
ical in the bulk and sphere-cap shaped for the hetero- 
geneous case, Vdrop becomes 4/37ri? 3 or 4/3irR 3 f(9), re- 
spectively, and we obtain R(p) and finally F S (R) and 
F s (R,0). 

Fig. S] a highlights the main results of our study. 
For the bulk we can directly test the classical theory 
of homogeneous nucleation by comparing our simula- 
tion results for F S (R) with the capillary approximation 
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FIG. 4: a) Plot of F 3 {R,9)/k B T vs. R for k B T/J = 3.0 and 
a broad range of values for 8. Dotted lines are results for 
L = D = 30, dashes lines for D = L = 40, dash-dotted lines 
for L = 60, D = 20. For the systems (bulk, ffi/J = and 
Hi/ J = 0.4 (6* « 58°)) the classical predictions 4tyR 2 j V £ and 
4-KR 2 -y v tf(9) are included (full curves), b) Plot of r vs. 8 
for fesT/J = 3.0. Crosses are data for D — L — 30, circles 
for D = L — 40, diamonds D — 10, L — 60, and squares for 
D = 10, R = 80. 



F c ,buik{R) = 47ri? 2 7^ 0, HI] (^ vl being the interfacial 
tension for a flat interface.) For R = 5 the difference be- 
tween the theoretical estimate and our simulation results 
is already less than 5% (at ksT/J = 3.0) and quickly 
vanishes with increasing droplet radius. The validity of 
the capillarity approximation for the bulk agrees with 
conclusions drawn from surface force measurements on 
liquid bridges [29(, but contrasts a large body of other 
work (see [26[ for a discussion). 

We have also determined the surface free energy for a 
sphere-cap shaped droplet attached to the wall for sev- 
eral contact angles as a function of the radius of cur- 
vature. It is important to note that within statistical 
errors different choices of L and D yield identical re- 
sults. For 6 = 90° and 9 = 58.2° we have also in- 
cluded the theoretical estimate from the capillarity ap- 
proximation F c j let = F c ,bulk ' /(#)■ When we consider 
Fc,het(R>@) ~ Fs(R,0), w e find that the difference in- 
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FIG. 5: Plot of the thermodynamic potential G for slab con- 
figurations (insert), normalized by temperature and total sur- 
face area 2LD, versus 2/D so that the slope yields the line 
tension (r/ksT = 0.26±0.001). The intercept of the ordinate 
was fixed to the literature value of the interface tension 0.434 

creases linearly with R, and hence can be interpreted 
as a line tension contribution, I-ett. I.e., Eq ^ is mod- 
ified to F s<het = F Syhom ■ f(9) + 2nrT. This division of 
surface and line contributions relies on our droplet defi- 
nition, Eq ([5]) , which uses the fact that for the lattice gas 
the Tolman length 3-5] is zero. Fig. HJd shows a plot of 
the line tension r obtained in this way versus contact an- 
gle. The line tension is negative, and becomes very small 
(and presumably vanishes) as 9 — > 0. The estimates for 
r have been used (together with the appropriate esti- 
mates for R) in Fig. [3J to show the expected deviation of 
the contact angle for small droplets from the bulk value. 
An independent check of the line tension is obtained for 
9 = 90° by investigating slab configurations (FigJSJ) with 
varying D. In this case the potential can be written as 
G/(2DLk B T) = ji v +2t/D, and the slope in Figf5]agrees 
with the estimate shown in FigHJ}. Hence, Eqs. (Q| and 
([5]) hold if a correction due to line tension is considered. 
Interestingly, an experimental hint for the need of a neg- 
ative line tension describing bubble formation in water 
extrusion from mesopores was given in [30| . 

In summary, we have shown that the classical the- 
ory of heterogeneous nucleation at flat walls, which pre- 
dicts a reduction of the free energy barrier by a fac- 
tor f(9), can describe the actual surface free energies 
F s (R, 9) , provided a line tension term is included into 
the latter. While the importance of line tension effects 
on small nuclei has been stressed in various other con- 
texts, e.g. 0, [I?], the present paper is the first ap- 
proach which provides a systematic method to obtain 
both F s (R) for spherical droplets in the bulk and F s (R, 9) 
for wall-attached droplets, as well as for the contact an- 
gle 9 and line tension r. Thus, for the first time we verify 
Turnbull's equation from 1950, which still stands as one 
of the cornerstones of the theory of heterogeneous nucle- 
ation. The methods, which in this letter are described 
for the simple Ising model, can readily be generalized to 



various models in broad classes of systems. They will 
enable significant progress in the understanding of nucle- 
ation phenomena in diverse branches of physics. 
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